function stars = countstars(c,s)

CV = abs(norminv([.10 .05 .01]./2));

stars = arrayfun(@(cv) abs(c./s)>=cv,CV,'unif',0);

stars = sum(cat(3,stars{:}),3);

stars = arrayfun(@(r) repmat('\ast',[1 r]),stars,'unif',0);

if numel(stars)==1
    stars = stars{1};
end